V us from 7r and K decay constants in full lattice QCD with physical u, d, s and c quarks 
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We determine the decay constants of the it and K mesons on gluon field configurations from 
the MILC collaboration including u, d, s and c quarks. We use three values of the lattice spac- 
ing and u/d quark masses going down to the physical value. We use the Wo parameter to fix 
the relative lattice spacing and f v to fix the overall scale. This allows us to obtain a value 
f° r /jf+//ir+ = 1-1916(21). Comparing to the ratio of experimental leptonic decay rates gives 
\Vus\ = 0.22564(28) Bl .( X +)(20)EM(40)iatt(5)v lld and the test of unitarity of the first row of the 
Cabibbo-Kobayashi-Maskawa matrix: \V ud \ 2 + \V US \ 2 + \V ub \ 2 - 1 = 0.00009(51). 



I. INTRODUCTION 

The annihilation of a charged tt or K meson to leptons 
via a W boson is a 'gold-plated' process whose rate can be 
determined very accurately from experiment. The decay 
width for pseudoscalar P made of valence quarks ab is 
given by: 

r(P i„) = °J^lf pm fM P (l " ||) 2 (1) 

up to known electromagnetic corrections. Here V a b is 
the appropriate Cabibbo-Kobayashi-Maskawa matrix el- 
ement for coupling to the W and fp is the pseudoscalar 
decay constant which parameterises the amplitude for the 
annihilation, fp can only be determined accurately from 
lattice QCD calculations. 

The experimental determination of T(K + — > 
Iv) /T(tt + — > Iv) can be converted to a result for the ratio 
of CKM element x decay constant for K and w. Using ex- 
perimental averages [JJ r(7r+ -> Iv) = 3.8408(7) x lO^" 1 
and T{K+ -> Iv) = 5.133(13) x lO^" 1 gives: 

= 0.27598(35) Br (A-+)(25) E M. (2) 

Here we have allowed for an electromagnetic correction 
to the ratio of widths given by (1 + <5em) with S^m — 
-0.0070(18)[1 The error in Eq. Q from this cor- 
rection is sizeable but not as large as that from the K + 
branching fraction to \xv which dominates. The total er- 
ror in determining the ratio of Eq. ([2| from experiment 
is then 0.16%. 

The electromagnetic correction means that fx and f w 
are defined as quantities in pure QCD without electro- 
magnetic interactions. An accurate theoretical result 
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from lattice QCD for f K + / f„+ then yields |Ks|/|Kd| U- 
Since V u d is known accurately from nuclear j3 decay, this 
gives V us . The higher the accuracy on V us the more strin- 
gent the test of CKM first row unitarity we can do, since 
V u b is too small to contribute (at present) to this. Any 
deviations are indications of new physics and the more 
stringent the test, the higher the scale to which the new 
physics is pushed. 

State-of-the-art lattice QCD calculations have 
achieved errors below 1% in fn/f-w [SHS], typically 
dominated by the systematic errors from extrapolation 
of the lattice results to the real-world continuum and 
chiral limits where the lattice spacing is zero and the 
u/d quark masses take their physical values (equivalent 
to the 7r meson mass taking its physical value). This 
means that significant improvements can be expected if 
we reduce discretisation errors, to make the continuum 
extrapolation more benign, and if we work with physical 
u/d quark masses that obviate the need for a chiral 
extrapolation. It also means that comparison of different 
methods for arriving at a physical answer from lattice 
QCD are important in testing systematic error estimates. 

The MILC collaboration recently gave an analysis 
°f ffc+/fir+ [10] from lattice QCD on their 'second- 
generation' gluon field configurations that include u, d, s 
and c quarks in the sea using the Highly Improved Stag- 
gered Quark (HISQ) formalism [TT] and a fully 0(a s a 2 ) 
improved gluon action I12j . They have ensembles with 
the average of the u and d quark masses down to the 
physical value. Their final error was 0.5% on the decay 
constant ratio, dominated by errors from the extrapola- 
tion to zero lattice spacing. Their analysis [TU] concen- 
trated on the ensembles with physical u/d quark mass 
and the aim was to perform a single self-contained anal- 
ysis that did not use additional information from, for 
example, chiral perturbation theory or determination of 
the lattice spacing using other quantities. 

In this paper we provide a new analysis with the most 
accurate result to date. To do this we use the same MILC 
ensembles with a completely independent analysis of me- 
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TABLE I: Details of the MILC gluon field ensembles used in 
this paper |14l 115] . /3 = 10/g 2 is the SU(3) gauge coupling 
and L and T give the length in the space and time directions 
for each lattice. ami :Sea ,am Staea and am CjSea are the light 
(up and down taken to have the same mass), strange and 
charm sea quark masses in lattice units. The ensembles 1, 
2 and 3 will be referred to in the text as "very coarse", 4, 
5 and 6 as "coarse" and 7 and 8 as "fine". The number of 
configurations that we have used in each ensemble is given in 
the final column. 



Set 


P 


ame tBe a 


QiTTLs ,sea 




L/a x T/a 




1 


5.80 


0.013 


0.065 


0.838 


16x48 


1020 


2 


5.80 


0.0064 


0.064 


0.828 


24x48 


1000 


3 


5.80 


0.00235 


0.0647 


0.831 


32x48 


1000 


4 


6.00 


0.0102 


0.0509 


0.635 


24x64 


1052 


5 


6.00 


0.00507 


0.0507 


0.628 


32x64 


1000 


6 


6.00 


0.00184 


0.0507 


0.628 


48x64 


1000 


7 


6.30 


0.0074 


0.0370 


0.440 


32x96 


1008 


8 


6.30 


0.0012 


0.0363 


0.432 


64x96 


621 



TABLE II: Values of the lattice spacing for the ensembles of 
Table [I] in units of parameters wo [13] . to [16] and r\ |17| . The 
ri/a values were calculated by MILC and given in [15] . 



Set 


w /a 


y/to/a 


ri/a 


1 


1.1119(10) 


1.0249(5) 


2.059(23) 


2 


1.1272(7) 


1.0319(3) 


2.073(13) 


3 


1.1367(5) 


1.0357(2) 


2.089(8) 


4 


1.3826(11) 


1.2389(5) 


2.575(17) 


5 


1.4029(9) 


1.2475(4) 


2.626(13) 


6 


1.4149(6) 


1.2521(3) 


2.608(8) 


7 


1.8869(39) 


1.6515(16) 


3.499(24) 


8 


1.9525(20) 


1.6769(7) 


3.565(13) 



son correlation functions with high statistics. We include 
ensembles with heavier-than-physical u/d quark masses 
and use chiral perturbation theory to pin down the point 
corresponding to physical light and strange quark masses. 
We also include very accurate information on the relative 
lattice spacings of the ensembles using the Wilson flow 
parameter, w [13] , This enables us to reduce the error 
on the decay constant ratio to below 0.2% which is close 
to the error coming from experiment in Eq. ([2]). 

Section [TTI describes the lattice calculation. The results 
and analysis, including a table of all our raw lattice values 
for meson masses and decay constants, are given in Sec- 
tion |III[ This is followed by a discussion and conclusions 
in Sections [IV] and [V] respectively. 



II. LATTICE CALCULATION 
A. Meson Correlators 

Table [I] gives the parameters for the MILC ensembles 
of gluon field configurations that we use here [21 ITS] . 
The table includes the values of sea quark masses in lat- 



tice units, where the u and d quarks are taken to have 
the same mass, mi = m u — m^. The accurate deter- 
mination of the lattice spacing will be discussed further 
below. Here we simply note that the 'very coarse' lat- 
tices (sets 1, 2 and 3) have lattice spacing, a « 0.15fm, 
the coarse lattices (sets 4, 5 and 6) have a ss 0.12 fm and 
the fine lattices (sets 7 and 8) have a «0.09 fm. Thus 
the spatial volumes of the lattices are large: the sets with 
physical mi (sets 3, 6 and 8) are all larger than 4.8 fm 
on a side, with sets 6 and 8 being larger than 5.5 fm on 
a side. 

On these ensembles we calculate light and s quark 
propagators using the same HISQ action as used in the 
sea. The valence t quarks are taken to have the same 
mass as those in the sea, the valence s quarks are re- 
tuned slightly to correspond more closely to the physical 
value [IB] ■ The valence masses used are given in Table III 
We use a unsmeared random wall source on each of 16 
time sources per configuration for very high statistical 
accuracy [15] , 

The propagators are combined to make meson corre- 
lation functions for 7r, K and rj s mesons. The rj s is a 
fictitious ss meson that is not allowed to decay here 
because we do not include the disconnected pieces of 
the correlation function. Since it does not contain va- 
lence u/d quarks it is a useful particle to study in lattice 
QCD [TU [19] . We include it here to provide more infor- 
mation to our fits about the meson mass dependence of 
the decay constants. 

We average over the time sources to obtain a result for 
each configuration and then study the configuration-to- 
configuration correlations. An autocorrelation analysis 
was performed in |18| and plots of the autocorrelation 
function for 7r and rj s correlators given for the ensembles 
with mi/m s = 0.1 and 0.2. The autocorrelation func- 
tion for the sets at physical mg/m s show a decrease in 
correlation between configurations. We bin over 2 adja- 
cent configurations on sets 1, 2, 3 and 6 and bin over 4 
adjacent configurations on all other sets. 

We then fit all three meson correlators simultaneously 
as a function of time, t, between source and sink accord- 
ing to: 



G mcson (t) = J^a k (e- E ^+e- E ^ T -^) (3) 

- (-i)*/-f;a fco ( e -^ t + e -^ T -*)). 

ko=0 

The oscillating piece is absent for 7r and rj s mesons be- 
cause the valence quark and antiquark have equal mass. 
We use Bayesian fitting methods (20] [21] so that the full 
effect of excitations in the spectrum can be included in 
the errors on the ground-state quantities that we are in- 
terested in, i.e. ao and Eq for each meson. The simulta- 
neous fit to all three mesons allows us to take into account 
the correlations between the fit results for each meson in 
our subsequent analysis. 
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(a/Wa) 



1 

0.98 
0.96 
0.94 
0.92 

0.9 
0.88 
0.86 
0.84 
0.82 

0.8 



ni|/m s =phys i — A- 



0.2 



0.4 



0.6 



0.8 



(a/w ) 



FIG. 1: At the top the ratio n /wo plotted as a function of a 2 
and below the ratio i/to/wo. The grey band gives the results 
of a simple polynomial fit to the a 2 and me dependence as 
described in the text. 



The 7r, K and rj s meson masses are given by the appro- 
priate Eq values from the fit above. The decay constant 
is determined from the corresponding amplitude, do, by 



fab = (m a + m b 




(4) 



for meson with quark content ab 5j. This formula holds 
for Goldstone pseudoscalar mesons made of staggered 
quarks and follows from the existence of a partially con- 
served axial current relation in these formalisms. The 
decay constant is then absolutely normalised in lattice 
QCD. 

Tabic |III| gives the results of our correlator fits for the 
decay constant and meson masses in lattice units on each 
ensemble. The errors are below 0.1% in almost all cases. 



B. Lattice Spacing Determination 

It has recently been proposed that lattice spacings 
be measured by smoothing the gluon field using a se- 
ries of infinitesmal 'smearing' steps [Ml [16]. This drives 
the gluon field towards a smooth renormalised field and 
gauge-invariant functions of this field, such as the pure 
gluon action density, become physical quantities. Param- 
eters used to determine the lattice spacing can then be 
defined from the flow-time, t, dependence of such quan- 
tities. Ref. pi)] defines to from: 



t 2 (E)\ 



t=t 



= 0.3 



(5) 



where (E) is the expectation value of the gluon action 
density. The parameter wq is preferred in [TJ], where Wq 
is defined by: 



0.3. 



(6) 



Wq should be less sensitive to small flow times where dis- 
cretisation effects may be important. Both wq and to can 
be determined by direct measurement on the gluon field 
and this makes them simpler to evaluate as well as typ- 
ically more precise than parameters based on the heavy 
quark potential [22 . The heavy quark potential must be 
determined by fitting large Wilson loops as a function 
of (lattice) time and then, to extract parameters such as 
r\ |14j . a further fit as a function of r must be done to 
the potential. 

None of wo, to and r\ can be simply related to any di- 
rectly measurable experimental quantity and their phys- 
ical value must be determined by a lattice QCD calcu- 
lation of such a quantity. For example, Wq and to are 
determined from the mass of the baryon in 13] and r\ 
is determined from a basket of quantities including the T 
excitation energy and the decay constant of the r\ s meson 
in [HUH]. 

Here we will use wo /a to determine the relative lattice 
spacing between the ensembles and finally fix its value 
from f n in our analysis in Sectio n |III| We give values for 
Wo/ a on each ensemble in Table [Il| and also, for compar- 
ison, of \/io/a. These were obtained using the methods 
explained in [T3J- We bin over 12 adjacent configurations 
to remove the effects of autocorrelations in the results. 

Fig. [T] compares the scales Wo, to and n by plotting 
r\/wo and y/to/wo as a function of (a/wo) 2 - We see 
that discretisation errors largely cancel between r\ and 
Wq since their ratio is very flat in a 2 . There is a small 
variation with mg. The errors here are dominated by 
statistical/fitting errors in r±. In contrast v^o/^o is 
much more precise but has relatively strong lattice spac- 
ing dependence, presumably from ^fio [T3]. We carried 
out a simple polynomial fit in me/ (10m s ) and a s (Aa) 2 , 
(Aa) 4 , (Aa) 6 with A=0.6 GeV for both ratios, taking pri- 
ors on the coefficients of 0.0(1.0). This gives the result 
n/w = 1.790(25) and y/to/w = 0.835(8) in the contin- 
uum and physical light quark mass limits, represented by 
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the grey shaded bands in Fig. [l] The result for y^o/^o 
agrees well with the BMW result of 0.835(15) (7) using 
the Wilson clover 2-HEX action in [IB] , 



III. ANALYSIS AND RESULTS 



Table |III| gives our raw lattice results for the pseu- 
doscalar meson masses and decay constants in lattice 
units. In this section we use these data to compute the 
dependence of /„., fx and f„ s on the quark masses and 
on the lattice spacing. This allows us to interpolate to 
the physical values of the strange-quark and light-quark 
masses, and to extrapolate to zero lattice spacing, ob- 
taining new predictions for the decay constants (as well 
as the rj s mass). The differences between our most chiral 
simulation data and our final results are small since our 
simulation is very close to physical. It is nevertheless im- 
portant to model these corrections accurately to optimise 
the precision of our final results. 

Our analysis involves the following steps: 

1. Remove the lattice spacing by multiplying the 
masses and decay constants for the it, K, and rj s 
by the values (with errors) of wo /a from Table [TTJ 
We also apply an svd cut to the data to guarantee 
that roundoff errors are not an issue when inverting 
the covariance matrix for the \ 2 function. This in 
effect triples the statistical errors. 

2. Fit the simulation results for wq/„, Wq/k, wofn B > 
and WqM^, together with the experimental result 
for /„.+ , as functions of the corresponding pion and 
kaon masses, and wq. We take the functional de- 
pendence from one-loop partially-quenched chiral 
perturbation theory plus terms polynomial in M%, 
M K , and a 2 . The fit gives a new value for wq , which 
is largely determined by the experimental value for 
f n used in the fit. It also gives the functional de- 
pendence of the decay constants and the rj s mass on 
the quark masses, as specified by M„ and Mjc ■ Sea 
and valence quark masses are specified separately. 
The same chiral formulas, with the same couplings, 
are used for pions, kaons and 77 s s; only the valence 
quark masses differ. 

3. Evaluate the best-fit functions for fx and f„ at val- 
ues of the pion and kaon masses appropriate for /„.+ 
and fx + . We set the u and d quark masses equal in 
our simulations. We correct for this approximation 
through appropriate choices for the values of M v 
and Mk used in our fit formulas to obtain our final 
results. At the same time we correct for (small) 
electromagnetic corrections to the meson masses. 
(The decay constants, by definition, do not need 
electromagnetic corrections; these are included ex- 
plicitly in Eq. Q.) 

In the rest of this section we elaborate on these steps, 
and survey our results. 



A. Chiral Fit 

We fit our lattice results for the decay constants using 
a formula drawn from partially-quenched chiral pertur- 
bation theory [23] that has the following form: 



/nLoK, X b , xl Ca , Xf a , L) + Sf x + Sf Ut . 



(7) 



Here /nlo is the result from chiral perturbation theory 
through 1-loop order, in a finite volume of size L on a 
side; and Sf x and <5/i a t are corrections for higher-order 
chiral contributions and nonzero lattice-spacing errors, 
respectively. We specify valence and sea quark masses 
through the dimensionless parameters x a , x^, etc. where, 
for example, a light-quark with mass nig = (m u + m^) /2 
would correspond to 



xi 



M % 
16tt2/o 2 ' 



(8) 



Here /o = 130 MeV is the experimental result for /„, 
and Mo )7r is the bare pion mass, obtained by subtracting 
the 1-loop chiral correction from masses measured in the 
simulation (Table III I . Using bare meson masses corrects 



for (negligible) finite-volume errors in the masses. The 
s-quark parameter is given by 



2M* 



ML 



(9) 



where Mq^k is the bare mass coming from the kaon 
masses measured in the simulation. The same formula 
is used for each of the three mesons we study, changing 
only the valence masses: 



u 

Sk 

u 



/nlo (a*, xi, xf a , xf a , L) + Sf x + 5f l&t 
hho(xe, x s ,x s e ca , xf a , L) + Sf x + dfu t 
f N Lo(x s ,x s ,x s r, xf a , L) + 6f x + Sfut (10) 



In our fits, we use very broad priors for the chiral pa- 
rameters in /nlo — 10-100 times wider than the final 
errors — so these have no impact on the fit. We also 
introduce a new parameter that multiplies the finite- 
volume correction in /nlo- This allows our fit to correct 
(crudely) for finite- volume corrections from higher orders 
in chiral perturbation theory. We set its prior to 1±0.33. 
Finite- volume corrections are quite small on almost all of 
the ensembles, as is evident from Tablc |TV] which lists cor- 
rections for the decay constants. (The rj s mass has very 
small corrections, similar in magnitude to those for fn B )- 
The finite- volume corrections agree with those calculated 
in |25j as well as with the finite- volume analysis in [10] . 

The square of the 77^ mass is fit with an analogous 
formula of the form 



M* LO (x a , x b , x s t, xT, L) + 8M* + SM, 



1 

lati 



(11) 



with x n 



5 



TABLE III: Values for it, K and rj a masses and decay constants in lattice units calculated for valence masses given in columns 
2 and 3. Some of the results were previously given in [181 123] . There are slight differences in some values with earlier results 
because different fit results were used. 



Set 




0,771 s ,val 






ciMk 


afK 


aM Vg 




1 


0.013 


0.0688 


0.23644(15) 


0.11184(10) 


0.41214(24) 


0.12695(15) 


0.53350(17) 


0.14185(9) 






0.0641 






0.40006(19) 


0.12585(10) 


0.51511(16) 


0.14009(7) 


2 


0.0064 


0.0679 


0.16614(7) 


0.10508(6) 


0.39077(10) 


0.12265(4) 


0.52798(9) 


0.14027(4) 






0.0636 






0.37948(10) 


0.12177(4) 


0.51080(9) 


0.13840(4) 


3 


0.00235 


0.0628 


0.10172(4) 


0.09938(6) 


0.36557(8) 


0.11837(4) 


0.50656(6) 


0.13720(2) 


4 


0.01044 


0.0522 


0.19158(9) 


0.09077(6) 


0.32789(11) 


0.10189(5) 


0.42358(11) 


0.11318(4) 


5 


0.00507 


0.0505 


0.13414(6) 


0.08452(5) 


0.30756(10) 


0.09788(4) 


0.41474(8) 


0.11119(3) 


6 


0.00184 


0.0507 


0.08154(2) 


0.07990(3) 


0.29843(5) 


0.09532(2) 


0.41478(4) 


0.11065(2) 


7 


0.0074 


0.0364 


0.14062(10) 


0.06618(5) 


0.23919(11) 


0.07424(4) 


0.30871(10) 


0.08236(3) 


8 


0.0012 


0.0360 


0.05716(2) 


0.05784(3) 


0.21855(5) 


0.06921(2) 


0.30480(4) 


0.08053(2) 



TABLE IV: Finite-volume corrections, A vo i/, to simulation results for the meson decay constants. Errors on the finite-volume 
correction come from our fit and are correlated between ensembles and between n and K. Also listed for each ensemble are the 
lattice spacing a (after determination of Wo which gives the error shown, correlated between ensembles), the ratio of valence 
strange to light-quark mass m 3 /me, the spatial dimension of the lattice L, and the pion and kaon masses (with their statistical 



a 


m s /me 


L 




M n 


M K 




A vo x/k: 




0.1543(8) fm 


5.3 


2.5 fm 


3.8 


302.4(2) MeV 


527.1(3) MeV 


1.24(23)% 


0.50(9)% 


0.10(0)% 


0.1522(8) fm 


10.6 


3.7 fm 


4.0 


215.5(1) MeV 


506.8(1) MeV 


0.38(7)% 


0.12(2)% 


0.00(0)% 


0.1509(8) fm 


26.7 


4.8 fm 


3.3 


133.0(1) MeV 


477.9(1) MeV 


0.43(8)% 


0.13(2)% 


0.00(0)% 


0.1241(7) fm 


5.0 


3.0 fm 


4.6 


304.5(1) MeV 


521.2(2) MeV 


0.37(7)% 


0.14(3)% 


0.01(0)% 


0.1223(6) fm 


10.0 


3.9 fm 


4.3 


216.5(1) MeV 


496.4(2) MeV 


0.24(5)% 


0.08(1)% 


0.00(0)% 


0.1212(6) fm 


27.6 


5.8 fm 


3.9 


132.7(0) MeV 


485.7(1) MeV 


0.15(3)% 


0.05(1)% 


0.00(0)% 


0.0907(5) fm 


4.9 


2.9 fm 


4.5 


306.1(2) MeV 


520.6(2) MeV 


0.41(8)% 


0.16(3)% 


0.02(0)% 


0.0879(4) fm 


30.0 


5.6 fm 


3.7 


128.4(0) MeV 


490.8(1) MeV 


0.21(4)% 


0.07(1)% 


0.00(0)% 



B. Higher-Order Corrections 

We include terms beyond one-loop order in chiral per- 
turbation theory by adding a correction of the form 

<5 f X = c 2a{x a + X b f + C 2b (x a ~ X h ) 2 

+ c 2c (x a + x b )(2x s r + xf a ) + c 2d (2x s r + xf a ) 2 
+ c 2e (2xf a2 + xf a2 ) 

+ C 3a (x a + Xbf + C 3b {x a + X b )(x a ~ X b ) 2 

+ c 3c (x a + x b ) 2 (2xT a + x s s ca ) 

+ Ci(x a + x b ) 4 + c 5 (x a + x b ) 5 + c 6 (x a + x b f 

(12) 

where we take priors of ± 1 for each parameter Cj . We 
only keep higher-order terms that might be significant 
given the precision of our simulation data. In fact, we 
obtain an excellent fit and almost identical results (to 
within a quarter of a standard deviation) when we keep 
only the quadratic terms. We include an analogous cor- 
rection, SM 2 , for the square of the r) s mass. 

We also correct for the nonzero lattice-spacing using 



Sfi 



lat 



aA 



QCD 



(13) 



where d n is allowed to depend upon the quark masses, 

d n = d n . + d nM {x a + x b ) + d n:lb (2xf a + x s s ca ) 

+ d nM (x a + x b ) 2 , (14) 

and again the coefficients have priors Oil. We get ex- 
cellent fits and almost identical answers without allowing 
for mass dependence, but we are conservative and include 
this possibility, since it could arise from taste-changing ef- 
fects [Tni US] , thereby increasing our final errors by about 
half a standard deviation. 



Eq. ( 13 ) is an expansion in the QCD scale Aqcd di- 
vided by the ultraviolet cutoff, n/a, for the lattice. The 
QCD scale is of order 500 MeV to 1 GeV. This is con- 
firmed by the Empirical Bayes criterion |20j which shows 
that our data imply a scale of about 600 MeV. In our 
analysis we use a conservative value, 1.8 GeV, to en- 
sure that nonzero lattice-spacing errors are not under- 
estimated. 



C. Isospin Violation and Electromagnetism 

We need to determine the xg and x s values correspond- 
ing to physical pion and kaon masses if we are to use our 
formulas to extract physical values for the decay con- 
stants and the rj s mass. The correct pion and kaon masses 



G 



come from experiment, but there are two complications 
that result from simplifications in the simulations. The 
first is that the simulation does not include electromag- 
netism. The second is that m u — rrid in the simulation, 
while in reality m u — 0.48(10) irid pQ. 

The most appropriate pion mass for /„.+ is the 
neutral-pion mass (134.9766(6) MeV [T]). All it mesons 
would have this mass in a world without electromag- 
netism — our simulations, for example — up to very small 
(quadratic) corrections from the u—d mass difference. 
These corrections are estimated at 0.32(20) MeV for M v + 
in [57] • For our purposes, it is sufficient to take 0.32 MeV 
as the uncertainty in the pion mass, and ignore the dis- 
tinction between charged and neutral pions: 



M ph ys = 134 .9 8 (32) MeV 



(15) 



This pion mass corresponds in our simulation to a 
light-quark mass of mi = (m u + rrid)/2- The corre- 
sponding kaon mass is one for an si meson. This is the 
root-mean-square average of the K + and K° masses with 
additional small corrections for electromagnetism: 



{Mf ys f 



M 2 K o) 



2 [(Ml + 

-(1 + A E )(M* + -M* )]. (16) 

Ae would be zero if electromagnetic effects in the K sys- 
tem mirrored those of the it. In fact it is closer to 1. 
Recent lattice calculations [28H30] that include electro- 
magnetic effects give values in the region 0.6-0.7. We 
take = 0.65(50) to conservatively encompass these 
results and this gives 



M 



phys 



494.6(3) MeV. 



(17) 



Tuning the pion mass to MP hys and the kaon mass to 
Mj^ hys in our fits sets the strange-quark mass to its phys- 
ical value, and the light-quark mass to the average mg of 
the u and d masses. This light-quark mass is correct, to 
within our errors, for the valence quarks in the pion, and 
for sea quarks in all three mesons. 

This tuning is not correct, however, for the 
K +, s valence light-quark, which is a u quark, with 
mass 0.65(9) rug. This difference produces a small but 
significant downward shift in fx+ ■ To compute the cor- 
rected K + decay constant, we evaluate our fit formulas 
with a pion mass given by -\/0.65(9)AfP hys , while adjust- 
ing the kaon mass so that 2M^ — m\ is unchanged (to 
leave the s-quark mass unchanged). These adjustments 
arc made only for the valence-quark masses in the K + ; 
the valence-quark masses in the pion and rj s , as speci- 
fied by M* s and Mjf ys , are left unchanged, as are the 
sea-quark masses in each of the mesons. 



D. Fit Results 



We fit wq times each of the decay constants and each 



0.1830 
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< 0.135 
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FIG. 2: Fit results for the it, K, and r/ s decay constants 
as functions of the light-quark mass for three different lat- 
tice spacings: 0.15fm (top/blue), 0.12 fm (middle/green), and 
0.09 fm (bottom/red). The data shown are from Table |Hl] 
with corrections for errors in the s masses, and for finite- 
volume errors. The lines show our fit with the best-fit values 
of the fit parameters. The dashed line is the a = extrap- 
olation, and the gray band shows our continuum results at 
the physical light quark mass point with nig — (m u + rrid)/2. 
The current experimental result for / T + is also shown (black 
point). Note that the three plots are against very differ- 
ent scales in the vertical direction: the range covered in the 
plot is 10 times larger than that covered in the f Vs plot. 



of the pion and kaon masses and wq. We also fit the 
experimental value for f n + — 130.4(2) MeV to our for- 
mula evaluated at the physical pion and kaon masses, 
17)). These fits are all done simultaneously 



Eqs. (15 



r] s mass in Table III to the formulas above, as functions 



using the same parameters for the fit functions in each 
case, and including the correlations between 7r, K and r/ s 
results discussed in Section HU 

The results for the decay constants, as a function of the 
light-quark mass, are shown in Figure [2j For each decay 
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TAB LE V: Sources of uncertainty in the final results 



0.000 0.005 0.010 0.015 0.020 0.025 

a 2 (fm 2 ) 



(Eq. ( 18 1) for the K + decay constant, the ratio of K + to tv + 
decay constants, the r/ s mass, and the Wilson flow parame- 
ter Wo- 
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chiral extrapolation 


0.03 


0.03 
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0.15 


a 2 — > extrapolation 


0.10 


0.10 
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0.27 


finite volume correction 


0.01 


0.01 


0.01 
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wq /a uncertainty 
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0.13 


0.03 


0.07 
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m-u/md uncertainty 


0.07 
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0.00 
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0.22% 


0.18% 


0.33% 


0.54% 



FIG. 3: Fit results for fic/fir evaluated at the physical light 
quark mass limit, with me = (m u +m [ j)/2, for different lattice 
spacings. The data shown are from Table [TTT] with corrections 
for errors in the quark masses, and for finite-volume errors. 
The top curve and data are from our analysis using wo to set 
the lattice spacing; the middle results are from our analysis 
using n instead of wo; and the bottom results are from our 
analysis using y/to. The gray band shows the final result from 
the wo analysis. 



constant we show fit results and simulation data for each 
of our three lattice spacings. The dashed line shows the 
continuum extrapolation, while the gray band shows our 
final results extrapolated to zero lattice spacing and the 
physical light quark mass limit (with the light-quark mass 
equal to the u—d average). The fit is excellent with a x 2 
per degree of freedom of 0.42, fitting 39 pieces of data. 
The final results are: 



U = 130.39(20) MeV 
f K + = 155.37(34) MeV 
f v , = 181.14(55) MeV 
M Vs = 688.5(2.2) MeV 

wo = 0.1715(9) fm 



/if + I /tt+ 

U/M Vs 
Ml/(2M 2 K -Ml) 



= 1.1916(21) 
= 0.2631(11) 
= 1.0063(64) 
= 0.9997(17) 

(18) 



Clearly the result for /„. contains no new information be- 
yond the input value from experiment that was included 
as a fit parameter. The K + results here are adjusted to 
correct the valence light-quark mass, as discussed above. 
We find that the K + decay constant is 0.27(7)% lower 
than the decay constant for a kaon whose valence light- 
quark's mass equals the u—d average mass. 

Error budgets for several of our results are presented 
in Table [V| Our fits are unchanged if we include addi- 
tional higher-order chiral or a 2 corrections, beyond what 
is discussed above. Omitting results from any one of our 
configuration sets shifts the mean values by no more than 
one standard deviation and usually much less. Omit- 
ting results from the smallest lattice spacing (0.09 fm) 
gives the same mean values but with standard deviations 
that are 2.5 times larger. Omitting the most chiral re- 
sults (m s /mi > 25) shifts the means by about 1/3 of a 
standard deviation and increases the standard deviation 



by 50%. These last two tests are evidence that our a 2 and 
chiral extrapolations are stable and robust. 

The a 2 variation of our simulation results is quite small 
(1-2 standard deviations) across our entire range of lat- 
tice spacings. This is illustrated in Fig. [3] where we show 
simulation results for fx/fn hi the physical light quark 
mass limit for our three lattice spacings (top curve); 
the gray band is the a = result. This behavior is in 
marked contrast with what we obtain if we set the lat- 
tice spacing using ^/Tq (bottom curve). The two methods 
agree to within 1.3 standard deviations when extrapo- 
lated to a = 0, but the variation with a 2 in the s/to anal- 
ysis is much larger. This agrees with the findings of [T3] 
that y/tg has larger discretisation errors than wq when 
compared to hadronic quantities. We have also redone 
our analysis using ri (middle curve). These results are 
similar to those from the wq analysis, and give an extrap- 
olated value that agrees with that analysis to within half 
a standard deviation. These two analyses also give: 



t = 0.1420(8) fm 
ri = 0.3112(30) fm (19) 

We use quite broad priors for wo, y/to, and r\ in our fits: 
0.1755(175), 0.1400(140), and 0.3150(320), respectively. 
They have little effect on the fit results. 

Finally, we give the Gasser-Leutwyler low-energy con- 
stants from the NLO term in our chiral fit. These are 
evaluated at scale M v and given in units of 10 -3 . 



L 4 = 0.36(34) 
L 5 = 2.00(25) 
2L 6 - Li = 0.28(17) 



L e = 0.32(20) 
L s = 0.77(15) 
2L 8 - L 5 = -0.46(20) 



(20) 

Values agree well with other chiral analyses, for example 
the MILC analysis on configurations including u, d and 
s asqtad sea quarks [31] . 



IV. DISCUSSION 



Eq. (181 lists a number of outputs from our analysis. 



The key result is that for f K + and in particular the ratio 
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15 1.17 1.19 1.21 1.23 1.25 



FIG. 4: A comparison of lattice QCD results for the ratio 
of K to 7r decay constants. The top two values (filled blue 
triangles) are for f K + //„.+ including u, d, s and c sea quarks 
- results from this paper and from [TO]. The lower values (red 
squares) include u, d and s sea and typically do not distinguish 
f K + from f K OT 



the continuum and physical light quark mass limits from 
those of the it and K . Here the surprising result, found 
earlier in [19] , is how closely the properties of the r] s 
match those expected from low order chiral perturbation 
theory. Our results here agree well with earlier results 
from n f = 2 + 1 [19] as well as from earlier analysis on 
these Tif = 2 + 1 + 1 configurations [IS] • 

We used wo/a to determine the relative lattice spac- 
ing between ensembles. Fixing the lattice spacing fi- 
nally from /„■ gives a physical value for wo in Eq. ( 18 1 of 
0.1715(9) fm. This agrees at 2a with the earlier result 
from BMW [13] of 0.1755(18) fm using the mass of the ft 
baryon to fix the physical value. The error in the BMW 
result is dominated by the statistical errors in the lattice 
calculation of Mq , and so we are able to obtain a smaller 
error using f^. 

In separate fits for comparison, not used for our ce ntra l 
values, we also obtained values for ^/to and r± in Eq. ( 19 1. 
Our y/to agrees with the value in [T3J. The value for 
r\ is not in good agreement with our earlier result on 
a subset of these ensembles [15] . however. The reason 
for this is largely because the values for r\ja have been 
updated and it underscores the difficulty of determining 
parameters from the heavy quark potential accurately. 



CONCLUSIONS 



Jk+ //tt+ needed to make use of Eq. This is obtained 
with an error of 0.18%. 

Fig. [4] compares our new result for fx+l f-n+ to earlier 
values on rif = 2 + 1 + 1 (u, d, s and c sea quarks) |10j 
and nf =2 + 1 (u, d, s sea quarks) configurations [SHS]. 
There is good agreement with earlier results within their 
larger error bands. Typically the results on nf = 2+1 did 
not distinguish between f K + and fx because the errors 
were not small enough to see this difference. The dif- 
ference that we see, 0.27(7)%, is in agreement with that 
expected from chiral perturbation theory (0.21(6)% [3]). 
It is also in agreement with a simple-minded argument 
assuming that the ratio of f Vs /fic depends linearly on 
Ml between and M% s . 

Our result for fx+l f-n+ is more accurate than MILC's 
recent analysis [TU] based on the physical point ensembles 
because we have included additional information: accu- 
rate relative lattice spacing values, fits to decay constants 
at heavier sea and valence quark masses, fits to 77 s masses 
and decay constants, and chiral perturbation theory to 
relate all of these fits at all of the lattice spacings to each 
other. The MILC error is dominated by their continuum 
extrapolation. Our Figs [2] and [3] show very little de- 
pendence on lattice spacing (when using wq) and benign 
extrapolations. Our error is nevertheless also dominated 
by the continuum extrapolation uncertainties along with 
statistical errors. 

Our analysis also gives results for the properties of 
the 77 s meson, by fixing its mass and decay constant in 



We give here the most accurate result to date for 
fx+ 1 f w + from lattice QCD. Our result comes from 
'second-generation' gluon field configurations with a 
highly improved discretisation of QCD and including u, 
d, s, and c quarks in the sea. We fit results from a range of 
ensembles with u/d quark masses down to the physical 
point and including additional accurate information on 
the relative lattice spacing between the ensembles. We 
test the robustness of our continuum extrapolation us- 
ing two different methods for lattice spacing determina- 
tion. These features mean that we are able to improve on 
the error obtained by the MILC collaboration 10 which 
aimed for a self-contained analysis using only physical 
u/d quark masses. 

Our result is 




With this level of accuracy the difference between fx 
with m u = rrid and fx + , which we determine to be 
0.27(7)%, is important. Going forward it will be nec- 
essary to make sure that f K and fx+ results from lattice 
QCD are averaged separately. 
Using Eq. ^ we determine 




= 0.23160(29) Br (K+)(21)EM(41) la tt. (22) 



It is no longer true that the lattice QCD error is much 
larger than the total of experiment plus corrections to 
experiment from electromagnet ism. 
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Given a value of V U d from nuclear j3 decay of 
0.97425(22) [32] gives 

\V US \ = 0.22564(28) Br(K+) (20) EM (40)i att (5)y w . (23) 

This agrees well with values from experimental results for 
semileptonic K decay rates combined with lattice QCD 
calculations of the appropriate hadronic form factor (2] 
121 [31] ■ The test of unitarity of the first row of the CKM 
matrix yields 1 - \V ud \ 2 - \V US \ 2 - \V ub \ 2 = -0.00009(51). 
This agrees well with the Standard Model result of zero 
and pushes the scale of new physics above 10 TeV [35] , 
To improve the limit on this scale significantly now needs 



further improvements to the accuracy of V u d |32j . 
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